Recent Advances in Dyson-Schwinger Studies* 



M.A. Pichowsky 
Center for Nuclear Research and Department of Physics 
Kent State University, Kent, OH 44242 U.S.A. 
(Dated: 11 February 2003) 

There have been many demonstrations of the utility of the Dyson-Schwinger equations of QCD 
as a systematic, phenomenological framework for describing the perturbative and non-perturbative 
dynamics of hadrons in terms of Euclidean Green functions of quarks and gluons. Still, there remain 
some unanswered questions regarding the theoretical underpinnings of the approach. I review several 
studies that are shedding light on how such questions involving the connection between the DSEs 
and their use in phenomenological applications might be resolved, and then review predictions for 
some exotic meson states. 
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I. INTRODUCTION 

The Dyson-Schwinger equations (DSEs) are an infi- 
nite set of coupled integral equations that relate all of 
the dressed n-point functions of a quantum field theory 
to each other. These n-point functions (which take their 
enumeration from the number of space-time points they 
depend upon) determine all of the dynamics of the theory. 
With no external particle sources present, the most ele- 
mentary of the n-point functions of QCD are the dressed 
quark and gluon propagators S(x\— x%) and D fll ,(xi— a^). 

A central aspect of QCD is that the elementary degrees 
of freedom appearing in the Lagrangian (quarks and glu- 
ons) are not simply related to the observed asymptotic 
degrees of freedom, the hadrons (mesons and baryons). 
This is referred to as confinement and understanding 
its origin represents one of the outstanding theoretical 
problems in physics. It then follows that knowledge 
of elementary two-point functions, the dressed quark 
and gluon propagators might not, by itself, be very 
useful since one can never observe quarks and gluons 
propagating directly. Clearly, n-point functions such as 
the fully-interacting quark-antiquark M(xi, . . . , X4) and 
three-quark scattering amplitudes H(xi, . . . , xq), which 
are related to meson and baryon observables, may be 
more interesting to explore. However, the self-coupled 
nature of the Dyson-Schwinger equations entails that 
these four-point and six-point functions are related to 
other n-point functions including the two-point dressed 
quark and gluon propagators S{x\— x%) and D^ v {x\—X2). 
Understanding the nature of this set of coupled, non- 
linear integral equations is a daunting task, but advances 
are continually being made. In the following, a summary 
of some of the most recent advances is provided. 

One traces the study of the DSEs of QCD back to 1950s 
and 1960's. These studies were crucial to understanding 
connections between field theories and observables, but 
applications to hadron phenomena were only moderately 
successful due to the lack of computational power avail- 
able at the time. In the 1980's, models formulated using 



the DSEs to describe the quark dynamics within hadrons, 
were found to provide good and compact descriptions of 
the light pseudoscalar and vector mesons. After these 
early phenomenological success, there followed a period 
of theoretical studies of the DSEs which explored the role 
of gluons, gauge invariance and confinement, which con- 
tinues to this day 0, Q ■ 

However, since the mid 1990's, much focus has turned 
to the study of physical processes which now permeates 
almost every subfield of hadron physics including: chiral 
physics and pion decays , 7T7t scattering electromag- 
netic form factors of the light pseudoscalar [5j and light- 
vector meson octets 0, inelastic scattering and quark- 
parton distributions |2j, high-energy diffractive electro- 
production and hadron scattering |8| . fin ite temperature 
studies of QCD Q, exotic mesons ^3, baryon bound 
states, their properties and electromagnetic form fac- 
tors > low-energy meson photo-production an d 
many others. 

Yet, a complete understanding of some theoretical as- 
pects of the framework has proven elusive. These unre- 
solved questions include: How does one construct sys- 
tematic and convergent expansions for the kernels of the 
DSEs? Are DSE and lattice QCD studies complemen- 
tary approaches to the study of QCD phenomena? How 
does one calculate observables in Minkowski space from 
a quantum field theory based in Euclidean space? In the 
following, I give a brief account of recent studies that pro- 
vide some answers to these questions. Finally, I report 
on an application of DSEs to exotic mesons. 



II. ROBUSTNESS OF TRUNCATION SCHEMES 



The simplest DSE is the Dyson equation, 
S" 1 ( P ) = Z 2 (C, A) i 1 ■ p + Z 4 (C, A) m(C) 



A j 4 

d q 



X a 
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~ Zl(C ' A) J j2^f g2D ^ ip ~ 9)y^s(g)r°(g,p),(i) 

which gives the dressed quark propagator S(p). In 
Eq. (UJ, A a for a = 1, . . . , 8 are the Gell-Mann SU(3)-color 
matrices, T"(q,p) is the dressed quark-gluon vertex and 
m(£) is the current quark mass. A translation-invariant 
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regularization scheme with scale A is employed to render 
the integral finite. The dressed quark propagator is of 
the form, 



S(p) 



Z(P 2 ) 



*7 ' P + M(j> 2 



(2) 



where Z(p 2 ) is the quark wave-function renormalization 
and M(p ) is the running quark mass. The dependence 
of Eq. on the ultraviolet scale A is removed by renor- 
malizing at p 2 — £ 2 , subject to the boundary conditions 
Z(( 2 ) = 1 and M(( 2 ) = m(Q, and absorbing the A- 
dependence into the multiplicative-renomalization con- 
stants Zi((, A). 

The dressed quark propagator is determined from 
Eq. only after the dressed gluon propagator D^ v (p—q) 
and dressed quark-gluon vertex T l v (q,p) are known. How- 
ever, these two- and three-point functions are solutions 
of their own DSEs, which in turn depend on still higher 
n-point functions. This is the infinitely-coupled nature 
of the DSEs, making it impossible to obtain a complete 
solution of QCD directly. 

In phenomenological applications, one may proceed by 
making assumptions regarding some subset of n-point 
functions such that the DSEs are reduced to a closed 
system of equations which may be solved directly. The 
simplest of such truncation schemes is the rainbow ap- 
proximation. It is equivalent to making the substitution, 



Z 1 (C,A)g 2 Tl(q,p) 



Pf 



(3) 



in Eq. QJ, where V(k 2 ) is an effective quark-gluon vertex. 
The exact dressed quark-gluon vertex depends on eight 
distinct Dirac-7 matrices and on three variables q 2 , p 2 
and q ■ p, while the effective vertex in Eq. maintains 
a single Dirac-7 matrix and depends on (q — p) 2 only. 

Assuming a given form for the dressed gluon propaga- 
tor D^(k) and effective vertex V(k 2 ), the quark propa- 
gator can be calculated from Eq. QJ. Then, one con- 
siders how a quark and antiquark combine to form a 
meson. This is described by the homogeneous Bethe- 
Salpeter equation (BSE), 

A(P 2 )r(p; P) = fS^Kip, q- P)S{q+)T(q- P)S(q-), 

(4) 

where K(q, p; P) is the two-particle irreducible qq scatter- 
ing kernel, T(j>; P) is the BS amplitude, and q± — q±P/2 
are the momenta of the quark and antiquark. This can 
be solved for any value of meson momentum P, but when 
the eigenvalue A(P 2 ) = 1, Eq. Q predicts the existence 
of a bound state meson with mass m = \J — P 2 . 

One can show that when the rainbow approximation 
is employed in Eq. Q), the ladder approximation for the 
Bcthc-Salpeter kernel, 

yi yi 

K(q,p;p) -» — 7 M (g) — 7 „ D^(q - p) V((q - p) 2 ), (5) 



must be employed in order to satisfy the axial-vector 
Ward-Takahashi identity and Goldstone theorem 0] . (A 
review may be found in Ref. Q-) Furthermore, it can be 
argued that the rainbow-ladder approximation represents 
the zeroth order term in a systematic series of approxi- 
mations that preserve the Goldstone theorem. However, 
what is truly amazing is that this series of truncations 
converges at a surprisingly rapid rate |13|, ! 



To illustrate how such truncations might converge, one 
compares solutions for the dressed quark propagator, 
and 7r- and p-meson masses obtained using the ladder- 
rainbow truncations to those obtained when higher-order 
terms are maintained in Eqs. © and J^J. In Ref. [T^j . this 
comparison was carried out using a simple model for the 
dressed gluon propagator D^ u (q — p). The quark-gluon 
vertex was dressed by a one-gluon loop, two-gluon loops 
and then with infmitely-many-gluon loops (using a recur- 
sive technique). The BSE kernel was constructed using 
the approach outlined in Ref. |l3[ . which ensured that 
the Goldstone theorem and axial- vector Ward-Takahashi 
identity were maintained at each level of the truncation. 



The study demonstrates that 7r and p meson bound 
state masses are well-described by the rainbow-ladder ap- 
proximation without significant change when dressing of 
the quark-gluon vertex is included [14|. (In retrospect, 
this explains the success enjoyed by the many studies 
of light pseudoscalar and vector mesons which have em- 
ployed the ladder-rainbow approximation.) The study 
also found that deeper in the time-like region, the dress- 
ing of the quark-gluon vertex may become more signifi- 
cant, but is still well- approximated by including only one- 
gluon or two-gluon loops; that is, results for the quark 
propagator and the BSE obtained by using the recur- 
sively dressed quark-gluon vertex were very similar to 
results obtained using just two-gluon loops to dress the 
quark-gluon vertex |l4| . 



The conclusion to be drawn from this study is that 
there is a systematic method for improving truncation 
schemes employed in phenomenological studies that al- 
lows a quantitative measure of their robustness. This 
achievement is quite important and represents a leap for- 
ward in the use of Dyson-Schwinger equations as a phe- 
nomenological tool. However, there is no free ride. Not 
all channels are equally robust or insensitive to the trun- 
cation schemes. Rather, pseudoscalar and vector chan- 
nels are particularly stable in this respect. There are 
other channels for which the inclusion of dressings aris- 
ing from a higher number of gluons has significant im- 
pact. Examples of such channels are the scalar and vector 
colored- diquark channels [Hil Il4|. Therefore, in practice 
one should check each channel of interest to verify that 
results are stable under changes to the level of truncation 
employed. 
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III. CONNECTION TO CALCULATION ON 
THE LATTICE 

The similarities between studies of DSEs and numeri- 
cal simulations of QCD on the lattice are clear [ll. ITol Il6| . 
Both explore QCD from its Lagrangian in terms of quarks 
and gluon degrees of freedom and both are formulated 
using the Euclidean metric where the square of the mo- 
mentum p 2 — p^p^ — p-p = Yli=i pf > is space-like. In 
principle, one might be able to compare results obtained 
from lattice QCD simulations to those obtained from the 
DSEs directly. In particular, lattice results may be used 
as input for DSE calculations or might provide further 
tests of the robustness of truncation schemes. Alterna- 
tively, Dyson-Schwinger studies may provide additional 
insights into the role of some symmetries damaged by 
lattice QCD simulations; for example, it is possible to 
use the Dyson equation (JJJ to extrapolate lattice results 
for the quark propagator to the chiral limit of vanishing 
quark mass. 

This was carried out recently |l7j | using a parametriza- 
tion of the Landau-gauge dressed gluon propagator 
Dfj, u (k) obtained from a lattice QCD simulation [l8j . and 
a simplified form of the dressed quark-gluon vertex of the 
form Eq. and parametrized by 

nt ^( 1 + ^P)( 1 + ^)-' + ..., m 

where logarithmic perturbative-QCD corrections have 
been suppressed, a and b are parameters determined 
by using this vertex and the lattice gluon propagator 
Dfj,v{k) in Eq. 0) to fit the dressed quark propaga- 
tor functions Z(p 2 ) and M{p 2 ). The resulting running 
quark mass functions M(p 2 ) are compared to those ob- 
tained from the lattice simulation [n| for a collection 
of current quark masses m(£) = 27, 50 and 102 MeV, 
as shown by the upper three curves in the left panel of 
Fig. 2] Self-consistency with lattice results required a 
dressed quark-gluon vertex that depends on the current 
quark mass, a(m) = A 2 (l + 7 (m/A g ) +125 (to/A,,) 2 )" 1 
and b = 0.005 A 2 , where £ = 19 GeV and A s = 
1020±100 MeV is a scale associated with the gluon prop- 
agator. 

The agreement between the mass functions calculated 
on the lattice and from the DSE JIJ is excellent for 
all finite values of the current quark mass to(£). Hav- 
ing parametrized the dependence on the current quark 
mass, we can then use this model to calculate the cor- 
responding dressed quark propagator in the chiral limit 
Z^, A)m(C) — > 0. The result is shown as the dotted 
curve in the left panel of Fig. Because of numerical 
difficulties, it is impossible at present for the lattice to 
simulate the chiral limit. However, a simple, linear ex- 
trapolation of the running quark mass function M(p 2 ) 
may be performed. The result of such an extrapola- 
tion is depicted by circles (lower data set) in Fig. ^ 
Clearly, the dotted curve fails to reproduce the linearly 
extrapolated lattice results for M(p 2 ). 
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FIG. 1: Comparisons of the running quark mass M(p 2 ) ob- 
tained from Eq. Q for current quark masses m(C) = 27 
(dot-dashed curve), 50 (dashed curve) and 102 MeV (solid 
curve) with corresponding results from lattice QCD simula- 
tions (data points) [Tg|. Also shown are M(p 2 ) in the chiral 
limit obtained from the Dyson equation (dotted curve) and 
from a linear extrapolation of the lattice results (circles). 

This disagreement is caused by having employed too 
simple a form for the dressed quark-gluon vertex in 
Eq. © , or from the inadequacy of a linear extrapolation 
of lattice data to the chiral limit. Using the Gell-Mann- 
Oakes-Renner relation as a check on self-consistency of 
the DSE approach suggests the latter may be true |17j . 
Reconciliation of these approaches requires additional in- 
vestigation, but it is clear that such comparisons are im- 
portant and will provide new insights for both Dyson- 
Schwinger and lattice communities in the future. 



IV. ANALYTIC CONTINUATION INTO 
MINKOWSKI SPACE 

Another issue, which has received less attention than 
those discussed above, concerns the mechanism by which 
a quantum field theory, formulated in Euclidean space, is 
related to hadron observables in Minkowski space. An- 
alytic continuation is the method by which observables 
in Minkowski space are obtained from the Euclidean am- 
plitudes or Green functions. For example, the p- meson 
bound state mass m p corresponds to finding the value of 
P 2 for which the BSE (2) eigenvalue \{P 2 = -to 2 ) = 1. 
Since the Dyson-Schwinger framework is defined for Eu- 
clidean momenta P 2 > 0, one must analytically con- 
tinue the Bethe-Salpeter eigenvalue A(P 2 ) to time- like 
momenta P 2 = —m 2 < 0. 

One way in which this process of analytic continuation 
may be understood is by allowing the fourth component 
of the Euclidean vector P M to become complex P\ — > iE 
where E > 0. As a result of this continuation, dur- 
ing the integration over the (Euclidean) relative-quark 
momentum q^, the arguments of the quark propagators 
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FIG. 2: The vector-meson ground state BSE eigenvalue 
A(P 2 ) for a range of meson momenta P 2 (solid curve). The 
meson mass m p is calculated using A(— m 2 p ) = 1, as indicated 
by the vertical line |2^|. 



q± = q 2 + iEqi — E 2 /4 sweep out a parabolic region of 
the complex momentum plane, whose maximum width 
is proportional to the meson's energy E — m p . To ac- 
complish this, one must know the behavior of the quark 
propagator S(q) within this parabolic region of the com- 
plex momentum plane. For light pseudoscalar (tt and K) 
mesons, the excursion into the complex plane is minor 
and so is easy to handle. For the light vector mesons 
(p, L), K* and <fi), this process can still be carried out 
by means of straight-forward techniques, albeit numeri- 
cally intensive |20j . However, for mesons with masses sig- 
nificantly greater than 1 GeV, one may encounter poles 
and/or cuts in the quark propagator during this process 
of analytic continuation. 

It has long been thought that the dressed quark prop- 
agator may contain complex-conjugate pairs of singular- 
ities and/or cuts in the complex-momentum plane |21| . 
To explore how one might carry out an analytic con- 
tinuation of the BSE in this event, we considered a 
simple model form for S(q) given by N pairs of complex- 
conjugate poles, 



S(q) 



N 

£ 

n=l 



«7 • q + m n ■ q + m* 



(7) 



where m n are complex-valued mass scales and z n are 
complex- valued magnitudes |22j. 

For the simple form given in Eq. and a phenomeno- 
logically successful BSE kernel 20], it was shown that 
one can deform the contour integrations thereby avoid- 
ing quark singularities, a necessary step to obtain the 
analytically continued eigenvalue A(P 2 ) for all values of 
momenta — oo < P 2 < +oo j^. The P 2 dependence 
of the resulting eigenvalue for the ground-state vector 
meson is shown as a solid curve in the right panel of 
Fig. El For comparison, the dashed curve in Fig. [3 was 
obtained by completely ignoring the need to deform the 
contour around poles when carrying out the integration 



in the BSE as this is not the correct implementa- 
tion of analytic continuation, the resulting dashed curve 
is not correct. The two curves are in perfect agreement 
for momenta P 2 > 4Re(m 2 ), where mi is the lightest 
quark mass parameter from Eq. JJJ. For such momenta, 
no quark singularities are encountered during the inte- 
gration and so, the curves coincide here. For momenta 
more time-like than this, the two curves diverge and the 
solid curve is the correct analytically-continued eigen- 
value A(P 2 ). The ground state vector meson mass m p 
is the solution of X(P 2 = —m 2 ) = 1. 

The Dyson-Schwinger framework is a renormalizable, 
Euclidean approach to the study of QCD. As such, one 
makes contact with observables by analytically contin- 
uing obtained amplitudes into Minkowski space. Using 
recently developed numerical methods, one can perform 
analytic continuations of the BSE, even in the presence of 
singularities in the quark propagators |22| . We now have 
the means to explore meson and baryon bound states of 
much heavier masses than previously possible. 



V. EXOTIC MESONS FROM THE 
BETHE-SALPETER EQUATION 

Before summarizing, I provide one example of a sim- 
ple model calculation of the meson spectrum. Studies of 
hadrons (both baryons and mesons) within the Dyson- 
Schwinger framework are too numerous to describe ade- 
quately in the space provided here. An extensive review, 
appearing recently, provides an excellent summary of the 
current status of such studies 3]. Rather than reiterate 
aspects of these studies, I give a brief accounting of a 
recent study of exotic mesons not covered therein. 

The most famous of exotic mesons is probably the 
7Ti(1400), the iso-vector, vector meson with positive 
charge-conjugation parity and quantum numbers J PC — 
1 '". The classification of "exotic" belongs to those 
mesons with quantum numbers that can not be con- 
structed as a quark-antiquark state . A quark-antiquark 
state \qq) with relative spin S and angular momentum 
L will have a spatial parity P = (— l) i+1 and charge- 
conjugation parity of C = (— l) i+s . To generate such a 
state, one must include additional degrees of freedom. In 
the quark model, this may be accomplished by introduc- 
ing either constituent gluons or flux-tube excitations or 
perhaps, by considering composite systems like \qqqq). 

In a quantum field theoretic framework, such particles 
are generated by constructing operators with the desired 
J = 1 h quantum numbers. One such operator is 
^i{x)^D^{x) where if)(x) is the quark field and is 
the fourth component of the covariant derivative = 
+ igA^ [10|. The first part of this operator may be 
used to project out the exotic meson bound state solution 
from the Bethe-Salpeter equation. Predictions for the 
masses of two exotic vector mesons were obtained |Toj 
using the simple separable model |23l ] which provides an 
excellent description of the light pseudoscalar, vector and 
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FIG. 3: Meson spectra obtained in Refs. IToll and I23I . One exotic pseudo-scalar with mass of 1319 MeV and two exotic vector 
mesons with masses of 1440 and 1499 MeV are obtained from the model. 



axial- vector meson spectra as shown in Fig. [21 One exotic 
pseudoscalar and two exotic vector mesons are predicted. 

VI. SUMMARY 

Recent theoretical advances are putting the Dyson- 
Schwinger framework on firmer ground and in closer 
contact with complementary frameworks, such as lattice 



QCD simulations. Future studies will continue to in- 
crease our ability to use this framework as a phenomeno- 
logical tool for exploring the dynamics of quarks and glu- 
ons in hadron processes. 
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